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Abstract 

In the framework of the nonrelativistic QCD and a quark-diquark model of baryons we 
have obtained the fragmentation functions for heavy quark to split into spin-1/2 and spin-3/2 
double heavy baryons. It was predicted the production rates as well as the shape of the energy 
spectra for the cc— and be— baryons in the region of Z° peak at LEP collider. 

as : 
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O ,. In the last decade the great success was obtained in the study of the heavy quarkonium production and decay (see, 
for example, Refs.Jl], |[). The investigation of the processes with the heavy quarks is based on the factorization 
hypothesis ||. The mass of heavy quark itiq is much larger than the scale of strong interactions Aqcd- So the 
presence of the small parameter Aqcd / m Q has allowed one to separate the effects of small and long distances. 
. The heavy hadron production amplitude may be presented as the product of the partonic part, which can be 
' calculated using the perturbative QCD, and the nonperturbative factor, which describes the free quarks to final 
ON . hadron transition. In the framework of the nonrelativistic quark model, this nonperturbative part can be presented 
through the quarkonium wave function at the origin \I/ (0) , which is calculated using the potential method (see 
, 0] and References therein). So, the approach based on PQCD and nonrelativistic quark model is known as a 
q ; nonrelativistic QCD (NR QCD) @. 

\Q . The heavy hadron production via fragmentation prevails in the large transverse momentum region at e + e~ and 

ON ' hadron colliders [||. Recently it was shown || that the fragmentation functions Dq^m(z, fi) and D g ^M(z, n) for 
quark and gluon to split into a heavy quarkonium with fraction z can be calculated within the framework of NR 
O .' QCD. The fragmentation functions are the process independent and can be applied to the e + e~, photonic and 
hadronic production of heavy quarkonia. Here we demonstrate that the NR QCD approach can be used for the 
calculation of the fragmentation function Dq^b(z, fi) for heavy quark to split into double heavy baryon . 

The estimation of the production rates for baryons containing two heavy quarks was done recently on the basis 
of the quark-hadron duality as well as on the basis of the PQCD ||. In spite of the some differences, it was 
suggested in Refs. that the double heavy baryon production has two step. In the first stage, there is the 
heavy quark Q fragmentation into double heavy diquark (QQ) in the colour octet state. The second step consists in 
the nonperturbative diquark fragmentation into a (QQq) baryon. It was mentioned in Ref. that this mechanism 
of the fragmentation can be factorized into short and long distance contributions as the same in Ref. §. The 
normalization of the fragmentation function Dq_>b(z,/z) is determined by the model dependent value of the (QQ) 
diquark wave function at the origin ^(qq)(0) as well as by the additional suggestion about nonperturbative diquark 
fragmentation mechanism into baryon. Note, that the calculation of the ^(0) for the colour object using potential 
model isn't grounded and destroys the factorization of the long and short distance effects. It was also suggested 
in Ref. that the probability of the diquark (QQ) fragmentation into baryon is equal to unity and the heavy 
diquark carries all of the momentum of the baryon -D(qq)^b(z, A 1 ) ~ ^(1 — z )- These circumstances raise too high 
the predicted values of the double heavy baryon production rates. 

The model of the quark-hadron duality, which was used in Ref. for the prediction of cc— baryon production 
cross section at B— factory, gives us the opportunity to obtain only upper limit of the production rates and doesn't 
predict baryon spectra. It was suggested in Q that all (cc)-pairs in the colour octet state with invariant mass 
from 2m c to 2Mq + AM (Mo is D-meson mass and AM = 0.5 -j- 1.0 GeV) split to the cc-baryon. However, it is 
known || that the phenomenological application of the quark-hadron duality method for heavy quarkonium photo- 
and hadroproduction needs to introduce additional process dependent if-factor which is equal to 1/3-^1/6 for 
description of the experimental data. Therefore the results of Refs. |?], ^] may be considered as a rough estimation 
of the double heavy baryon production rates, which need in more detail analysis. 
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2 The model 



In the present paper we consider a new mechanism of the double heavy baryon production via heavy quark 
fragmentation based on the hypothesis of the point-like diq uarks, which can be produced directly in the hard 
interactions of quarks and gluons at the short distance jnj. As it was shown in Ref. the data on e + e~ 

annihilation into hadrons don't contradict to the existence of the very small diquark (Qq) containing one heavy 
and one light quark. In such a way, the heavy quark (Q = b, c) can fragment directly into double heavy baryon 
spin-1/2 or -3/2 catching the scalar or the vector heavy diquark (D = (eg), q — u,d,s), correspondingly. In the 
process of the production quark-diquark system Q(Qq) in the colour singlet state the relevant QCD scale satisfies 
to the next condition: /i > ttiq + 2m o S> Aqcd, where mo is the heavy diquark mass. This fact leads to the 
factorization of the double heavy baryon production amplitude. The transition of the Q(Qq) into final baryon can 
be described using the nonrelativistic approximation of the potential model, because the system Q(Qq) contains 
two heavy particles. 

So, for the calculation of the baryon (B = Q(Qq)) wave function at the origin we have used equation p2| , which 
effectively takes into account the relativistic effects of the kinematic nature p3[: 



+ U QD (r))^{f)=E^(r), (1) 



V 2 
where 

p 2 M 2_ (TO |_ m ^2 

£= W 4AP ' 

p2 = [M 2 - (mq + m D ) 2 ][M 2 - (m Q - m D ) 2 ] 
AM 2 

The potential of the quark-diquark interaction is taken in the conventional form 0] : 

b 

Uqd = — + or, 2 
r 

where a — 0.183 GeV 2 , b — 0.52. The results of our numerical calculations for the radial part of the wave function 
at the origin |i?(0)| 2 = 47r|*(0)| 2 , are presented in Tables 1 and 2. 

The fragmentation function Dq^b(z, fi) at the scale \x = \x Q = niQ + 2m]j for the production of the baryon B, 
containing the heavy quark Q and the heavy diquark (Qq), is given by the next expression: 



1 f°° \M\ 2 
Dq^b[Z,h ) = / ds lim ==, (3) 

16?H J Sm . n 9o^oo |_/v^ |2 

where M. is the matrix element for the production of a baryon B and antidiquark D with the total four-momentum 
q and the invariant mass s = q 2 , Ai is the matrix element for the production of a quark Q with the same 
three-momentum q. The lower limit in the integral (3) is 

_ M 2 +f T _ ml+ft 
z 1 — z 

Here M — vtlq +mj is the baryon mass, pr is the baryon transverse momentum in the reference frame where 

Po + P3 



q = (g o ,0,0,g 3 ), p = (po,PT,P3) and z 
In the limit of q Q — > oo: 



q a + g.3 



_ Mr 



Z 1 — Z 

In the axial gauge associated with four- vector n = (1, 0, — 1): 



(kn) ' 



the fragmentation contribution comes only from the Feynman diagram for Q — > QDD shown in Fig.l. 
The gluon couplings to scalar and vector diquarks are presented by the next expressions: 

Sl = -ig s T\q' -p^^Fsik 2 ), (5) 
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Vl = ig s T b \ ehehiq'-PD)^ 2 ) - [tfe* D )e* Dlt - (p D e* D )s* Dtl }F 2 (k 2 ) 



(eU'WdPdW -p D )„F 3 (k 2 ) 



(6) 



where T b are Gell-Mann matrices, e* Dl e*j are the diquark polarization vectors. F S ,F\,F2 and F3 are form factors 
depending on the momentum transfer squared k 2 = (q' + pjj) 2 . 

The light diquark form factors may be parametrized as follows |l5[ : 



F S (Q 2 ) 



Ql 



Ql + Q 2 



Fi(Q 2 ) = 



Ql 
Ql + Q 1 



(7) 



F 2 (Q 2 ) = (1 + k)F 1 (Q 2 ), F 3 (Q 2 ) = 

where k ~ 1.39 being the anomalous magnetic moment of the vector diquark, Q 2 w Q 2 /2 and Q 2 « 3 GeV 2 Era. 

At present time the form factors of the heavy diquarks are unknown and one has to assume a certain dependence 
on Q 2 . The general consideration (see, for example, Ref. 11), based on asymptotic QCD picture and the 
phenomenology, predicts F s ~ 1/Q 2 , F% ~ Fj, Q 2 ^> 3 GeV^, corresponding to a very small heavy diquark. 
We try to take these points into account by using the simplest possible expressions for the heavy diquark form 
factors: 



F S (Q 2 ) = gf, 



F 2 (Q 2 ) = (1 + K )F 1 (Q 2 ), F 3 (Q 2 ) = 0, 



(8) 



where Ql = fc 2 air = 4m 2 ;,, Q 2 = Q 2 /2. It has k^ in « 20 GeV 2 for (cq) diquarks and our parametrization is the 
same as in Ref. (TTJ . Note, that such choice of the diquark form factors gives opportunity to obtain the compact 
analytical expressions for the fragmentation functions. 

In the case of the heavy quark fragmentation into spin-1/2 baryon it has fusion of the heavy quark Q and scalar 
diquark D. After the some obvious simplifications we have obtained: 



Mi /2 = *(0) 1 



M 



2m Q m D 3^3 (s • 



. A5 ij F s (k 2 ) 
72^2 



(9) 



2U( PQ ) 



-M{q + m Q ) + (s-m 2 Q ) 



(np) 
ink) 



where g s = \J 47ra s , pq — fp and po = °rp are the momenta of the quark and the diquark in the baryon, corre- 
spondingly, r — mu/M, f = 1 — r, p is the baryon momentum, 4<5 i:, /3"\/3 is the colour factor of the amplitude 
shown in Fig.l, the Dirac spinor T is the matrix element for the production of a Q quark of momentum q = p + q' , 
q' is the antidiquark momentum, the spinor U(p) describes spin-1/2 baryon in the final state. 

The amplitude for heavy quark fragmentation into spin-3/2 baryon, corresponding to fusion of the heavy quark 
and the vector diquark, can be written as a sum of two parts, which are proportional to Fi and F% form factors: 



M 



3/2 



M\ /2 



■^3/2. 



(10) 



M\ /2 = *(0) 1 



M 



.46*' F x (k 2 ) 
2~\2 



2mQm D s 3V3 (s - m 2 Q 



-M(q + m Q ) + (s - m Q ) 



2 \ (np) 



ink) 



M 



2mQm D s 3^3 (s - m 2 Q 



(11) 



A5 ij F 2 (k 2 ) 

2 \2 



(12) 
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D\ 



(q + m Q )T. 



Here, the spin- vector describes the spin-3/2 baryon in the final state and satisfies to the following conditions 

(p-Jlf)¥ M (p) = 0, 
= 2M, 7^ = P/.*/* = 0- 

In the case of the unpolarized spin-3/2 baryon the summing on the helicity states is carried out by means of the 
following formula: 

2^ *£(P)*£(P) = (P+ M ) [9^ ~ 3 7^ 



3M 2 



3M 



(13) 



3 The results 

Substituting the amplitudes (9) and (10) to the basic formula (3) we have obtained the fragmentation functions 
for a heavy quark split into the spin-1/2 and spin-3/2 baryons. Omitting the details of the calculation we write 
here final result for the heavy quark fragmentation function for spin-1/2 baryon: 



_. , . 8a 2 s (2m D ) |*(0)| 2 Q* , 

1/2 V >H-oJ 4Q5r 3 M 3 M4 1/2V , 



where 



z 4 (l-z) 3 



15 - 6z(7 + 3r) + z 2 (44 + 6r + 15r 2 ) 



-22z 3 (l-r)+5z 4 (l-r) 2 

The fragmentation probability for the production of the spin-1/2 baryon is 

f 1 „ / * . 8a 2 (2m D ) |f(0)| 2 Q 4 . . 
/ D- l / 2 (z,a )dz — — — — ^— L1 — ^-^--^Ii /o(r), 
J 1/2V ' 405r 3 M 3 M 4 1/2V ; ' 

1 



where 



h/i{r) 



\l - 12r + 75r 2 - 420r 3 - 1827r 4 - 126r 5 + 



42r 4 (l - r) 10 
2037r 6 + 300r 7 - 30r 8 + 2r 9 - 1260r 4 (l + r) 2 ln(r) 
The fragmentation function for Q split into spin-3/2 baryon is 

^3/2^,/io) = 1215r r M 3 MS^^r), 



where 



$n(z,r) 



F 3 /2(z, r) = + 4(1 + «)$ 12 + (1 + k) 2 $ 22 , 

z 4 (l-z) 3 



7z 8 (5r b - 18r 5 + 51r 4 - 76r 3 + 51r 2 - 18r + 5) 



(1-z + rz) 141 

f 2z 7 (105r 5 - 455r 4 + 1122r 3 - 1038r 2 + 413r - 147) + z 6 (105r 6 - 42r 5 + 1596r 4 - 
2648r 3 + 5581r 2 - 2198r + 1134) + 2z 5 (-21r 5 - 945r 4 - 100r 3 - 4352r 2 + 1505r 
-1295) + z 4 (847r 4 + 2028r 3 + 8031r 2 - 2170r + 3780) + 2z 3 (-446r 3 - 
2018r 2 + 343r - 1785) + 7z 2 (121r 2 + 2r + 302) + 42z(-r - 17) + 105], 



$22 (z,r) 



z 4 (l~z) 3 
(1 — z + rz) 



in 



[45 - 18z(7 + 3r) + z 2 (164 + 84r + 45r 2 



(14) 



(15) 



(16) 
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2z 3 (65 + 43r + 15r 2 ) + z 4 (47 + 56r + 65r 2 )], 
*ia(*,r) = ~ Z ^ 12 [z 6 (l - 9r + 17r 2 - 19r 3 + 10r 4 ) 




z 5 (7 - 28r + 35r 2 - 22r 3 ) + 3z 4 (ll + r + 22r 2 + 16r 3 + 10r 4 ) 
z 3 (82 + 87r + 95r 2 + 51r 3 ) + z 2 (103 + 98r + 47r 2 ) 
-3z(21 + llr) + 15]. 



The fragmentation probability for the production of the spin-3/2 baryon can be calculated in the same way as 
for spin-1/2 baryon. We don't present this result here because of it's unwieldy. 

The fragmentation function Dq^b(z, /u) satisfies to the Gribov-Lipatov-Altarelli-Parisi (GLAP) evolution equa- 



The boundary condition on the evolution equation is the initial fragmentation function Z)q_>b(z, /i c ) at the scale 
fi Q = rriQ + 2r7j£). Note, that at leading order in a s one has: 



and the evolution equation implies that the fragmentation probability 

/ Dq_>b{z, u 2 ) does not evolve with the scale /z. Therefore the fragmentation probability is the universal charac- 
teristic of the production rates. The evolution only changes the z-distribution to smaller values of z. 

The results of calculation for the fragmentation probabilities and the average values of the momentum fraction 
for Q quark to split into double heavy spin-1/2 and spin-3/2 baryons are shown in Tables 1 and 2. The values 
of < z > are presented at fj, = fi a and /j, = Mz/2. We used the following set of the mass parameters: m c = 1.7 
GeV, m& = 5.1 GeV, mf u = mf d = 1.9 GeV, mf s = 2.0 GeV. The masses of the vector diquarks are greater than 
corresponding scalar diquarks by 0.1 GeV. The fragmentation functions £> c ^e cc (z, /z) and Db^A bc [z, /i), divided 
by the fragmentation probabilities, are shown in Fig. 2 at fi = ji and fj, = Mz/2. In the range of accuracy of our 
model the shape of z— spectra for spin-1/2 and spin-3/2 baryons are the same one to other as well as z— spectra 
for baryons with strange quark (fl cc ,Qbc) and without it (S CC ,A(, C ). 

Let compare our results with the estimation [|| , where it was founded that fragmentation probabilities for heavy 
quark to split into double heavy baryons are (2 -=-3) • 10 -5 independently on the baryon spin and the flavour content. 
In our approach the fragmentation probability depends on baryon spin and flavour. 

Moreover, there is the relative growth of the production rates of the spin-3/2 baryons for 6-quark fragmentation 
in comparison with the fragmentation of c-quark. We have obtained that S* c /S cc ~ 0.7 and Aj* c /A;, c ~ 1.8. Our 
results for the fragmentation probabilities for the spin-3/2 baryons are approximately equal to the results of 
however the fragmentation probabilities for spin-1/2 baryons are smaller about factor 2. Taking into consideration 
the uncertainties of our calculation, connected with the diquark form factors (8) , we can predict the fragmentation 
probabilities for the cc— and be— baryons are about of 10~ 5 . There are more accurate predictions for baryon z- 
spectra in the our approach, which are practically independent on the diquark form factors and other parameters. 
We have obtained that the average fraction of the baryon momentum < z 0.54 for the cc— baryons and 
< z >« 0.66 for the be— baryons at fi = Mz/2, and < z >ps 0.75 and < z >ps 0.84, correspondingly, at [i = fi a . 
This results are independent from the baryon spin. By contrast, in the model, based on Q — > (QQ) fragmentation 
mechanism corresponding values at fi = [i Q are equal: < z >= 0.57-^0.62 for the cc— baryons and < z >= 0.68-f-0.73 
for the be— baryons. 

In conclusion, we sum the results of our paper. In the framework of NR QCD and a quark-diquark model of 
baryons we have obtained in the leading order in a s the fragmentation functions and the probabilities for heavy 
quark to split into double heavy spin-1/2 and spin-3/2 baryons. Using the QCD evolution equations we have 
recalculated the fragmentation functions from the initial scale fi to fi — Mz/2. Our results can be used for 
prediction of the double heavy baryon production rates as well as for the description of the energy spectra at 
e + e~-collider LEP. Contrary to Refs. |0 | we have predicted also few precise effects: the nontrivial spin and 
flavour dependence of the baryon production rates, the value of the average baryon momentum < z > and it's 
dependence from the scale /j, as well as from the baryon flavours. 



tions 




(17) 



where 
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Table 1. The spin- 1/2 baryons. 



Bar yon 


M, 
[GcV] 


\R(0)f, 
[GeV 3 ] 


Pq^b 


< z > 


< Z > Mz /2 


1 — 'CC 


3.6 


1.2 


5.8 • 10~ b 


0.75 


0.55 




3.7 


1.3 


5.3 • 10-" 


0.74 


0.54 


Abe 


7.0 


3.8 


8.9 • 10-" 


0.84 


0.67 


^bc 


7.1 


4.0 


7.4- 10~ 6 


0.83 


0.66 
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Table 2. The spin-3/2 baryons. 



Baryon 


M, 
[GeV] 


|iJ(0)| a , 
[GeV 3 ] 


Pq^b 


< Z >o 


< Z >M z /2 




3.7 


1.25 


1.2 • 10~ 5 


0.75 


0.54 




3.8 


1.35 


l.i • icr 5 


0.74 


0.53 


A* 


7.1 


4.0 


3.5 • 10" b 


0.84 


0.67 


o* 


7.2 


4.1 


3.3 • l(T b 


0.83 


0.66 



Figure captions. 

1 . Diagram used for description of the heavy quark to split into double heavy baryon. 

2. The fragmentation functions normalized to the unity at n = ji a (curves 1 and 3) and /i = Mz/2 (curves 2 
and 4). The curves 1 and 2 - Ds cc (z,fM), the curves 3 and 4 - D\ bc (z, fi). 
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